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Abstract. A non-iterative topological sensitivity framework for guaranteed far field detection 
of a dielectric inclusion is presented. The cases of single and multiple measurements of 
the electric far field scattering amplitude at a fixed frequency are taken into account. The 
performance of the algorithm is analyzed theoretically in terms of resolution, stability, and 
signal-to-noise ratio. 


1. Introduction 

A thriving interest has been shown in topological sensitivity frameworks to procure solutions of 
assorted inverse problems especially for detecting small inhomogeneities and cracks embedded 
in homogeneous media [a El 0 HU], The impetus behind this curiosity is the robustness 
and simplicity of these algorithms. However, the understanding and utilization was somehow 
heuristic. For the first time, Ammari et al. rigorously debated these algorithms [HEiEi]. 

Our aim is to design and debate a far field detection algorithm based on topological sensitivity 
of small dielectric inclusion Dp := pOo C with position z^i, scale factor p and a smooth 
bounded domain Od C The inclusion Dp (with permittivity ei = e^eir > 0 and permeability 
po > 0 ) is assumed to be embedded in (with background permittivity cq > 0 and permeability 
po > 0, setting forth a non-permeable behavior (but the results extend to the case otherwise). 
Here ei^ > 0 is the relative permittivity. The medium is probed by the incident electric field 


Eo(x) := Vx X = inO x x G (1) 

where k = oj^/IqPq G M+ is the wavenumber, a; > 0 is the angular frequency, 0 G := {x G 
I X • X = 1} and 9^ is any vector orthogonal to 9. Then, the total electric field Ep in in 
the presence of Dp satisfies 

V X V X Ep(x) - K^(l - (1 - ei^)XD^(x))Ep(x) = 0, x G (2) 

lim |x| [V X (Ep — Eq) X X — zk (Ep — Eq)] = 0, (3) 

|x|—>-oo 



where XDpi'^) is the characteristic function of Dp and v := v/|v| for all v G \ {0}. Let 
:= {v G L^(S^) | v(x) • x = 0, Vx G S^} (called far-field amplitude) be defined by 

Ep(x) - Eo(x) = ^x; + O (l/|x|^) , as |x|oo. (4) 

Then, the following problems are dealt with henceforth. 

(i) Given the single data set {E^(x; 0), Vx G S^, 9 G S^}, find location zd of inclusion Dp. 

(ii) Given the multiple data set {E“(x; 9j), Vx G §^, (0j)j=i G n G N}, find location Z£). 

In Section [2l a topological sensitivity based location search algorithm is presented and its 
resolution is ascertained. A stability analysis is performed in Section [3l 


2. Topological sensitivity framework and resolntion analysis 

Consider a search point zs G and create a trial inclusion Ds = zs + 50s inside with 
permittivity 62 = eoe 2 r > 0, scale factor <5 > 0 and smooth domain Os C Let E 5 be the 

electric held in in the presence of Dg subject to incidence Eg and satisfying ©-([S]) with 
{^ir,XDp, Dp) replaced with {e 2 r,XDs^ Ds). Let E^ be the far held amplitude of E 5 . 

Dehne the misht j7’[Eo] and its topological derivative (imaging) functional by 


J[Eol(»s) := f / 

4 J§2 


E- 


(x;0,0^ 


Er(x;0,0^) 


2 

ds(x), 


(9tJ’[Eo](z5) := 


dJ[Eo](z5) ,,, 
5((f)3 • ^ ’ 


It is claimed that zs G R^ relative to which dT,T[Eo] has a sharp peak is a candidate for zd. 
The claim is justihed in the sequel. Recall hrst the asymptotic expansion of E^ versus p [5]: 

E“(k; 9) = -K^p^iATT)-^ (er; - l) (l - MpEo(zz?)e-'""^"« + 0{p^), (6) 

wherein Mp = Mp(ej“^^,O d) is the polarization tensor. Note that Mp = 3eir(eir + 2)“^|Oz)|I 
for a spherical domain, where I is the 3x3 identity matrix and superposed T indicates matrix 
transpose. It is worthwhile precising that E“ also admits an analogous expansion versus 5. 
Further, consider the dyadic Green’s function 

r(x,y) := -eo [l + e'"l"-yl/47r|x - y|, (7) 

satisfying Vx x Vx x r(x, y) — K;^r(x, y) = —eo5y(x)I for all x, y G R^ subject to the Silver- 
Miiller condition. Here 5y(-) = <5o(' ~ y) is the Dirac mass at y and Vx acts on matrices 
column-wise. Moreover, for isotropic materials (see, for example, a) 


r(x,y) = r(y,x) and Vx x r(x, y) = [Vy x r(y, x)]'^ , Vx,yGR^. ( 8 ) 

Finally, recall that the electric Herglotz wave for all G L|(S^) is dehned by ([HI Ch. 6 ]) 


iRi?[$](z) := / $(x)e*'"^^^ds(k), z G R^. 
Js2 


(9) 


Lemma 2.1. For all search points zg G R^ 


5r77[Eo](z5) = -k^{At:)-^ - l) 3?e {iR£;[E-(-; 0)](zs) • M 5 Eo(zs)} 


( 10 ) 






Proof. Expanding misfit ([5|) and invoking expansion ([6]) for render 


J[Eo]{zs) - |E-(x;0)|'ds(x) = -3fte E^(x; 0) • E-(x; 0)ds(x)| + 

- l) 5Re| J^ (l - xx^) M 5 Eo(zd) • E“(x; e)e~ 






ds{k)\+0{6^). ( 11 ) 


Remark that for any matrix A, vectors u, v, density ‘h G Lf{E>^) and x, z G 


Au • V = A^v • u, 

X (x X <h(x)) e" 


(l — xx"^) u = —X X (x X u) 


= -Vz X Vz X <h(x)e' 


(ch(l 


( 12 ) 

(13) 


Therefore, by the fact that V x V x —k'^T-Le[^] = 0 for <h G Tf (S^) and by above identities 


J[Eo]{zs)-\ljB^{x-,e)\^ds{ 


2^3 


X = 


dvr 


€2^ - l)5Re {iRe[E-](z5) • M5Eo(z5)} + 0 ( 6 ^). 


Finally, the proof is completed by taking limit as 5^ ^ 0. □ 

In order to substantiate the capability of 9 t^[Eo] to locate Z£), we elaborate ^e[E“](z) as 

nE[E^i-; 9 )]iz)= [ E-(x; 0 )e--"^ds(x), 


/S 2 


^ ^ -l) [ (I - xx'^) ^°'>MpEo{zD)ds{x) + 0{p‘^), 

Js 2 


dvr 


k'^o^ 


1 T 

I--VzV^ 

K 




- l) 9m|r(z,Z£i)|MpEo(z£,) + 0 (p^). 

Indeed, the last identity is evident since 

f e*'^^^*-^“^°^(is(x) = jo(«:|z — ze\) = d7rK“^9'm{e*'^l^~^°l/47r|z — ze\], 
J §2 

where jn is the order n spherical Bessel function of first kind and consequently. 


MpEo(zD) + 0 (p^) 
(Id) 


I-^VzVj 


^ f =-—$5m|r(z,ZD)|. 

/ Js2 KeQ 1 J 


(15) 


(16) 


In fact, the following result is proved by virtue of Lemma l2.II and expansion (jldli . 

Theorem 2.2. For all zs G and constant ■= (e^^^ — l) — l) 

dTJ[Eo]{zs) = V^"C’.(4vreo)-'5Re{9m{r(zs,z,5)}MpE;)(i^-M5Eo(zs)} + 0(p^). (17) 
Remark 2.3. It immediately follows from Theorem \2.S\ that 


5 t 77 [Eo](z 5 ) oc 9m{r(z5,ZD)} = - 


dvr 


\jo{Kr)l + j2{nr) ( rr'^ - ^1 
















where r := zs — zd with r := |r| . Since jnikr) = 0{l/kr) as kr ^ oo and jnikr) = 0{{kr)'^) as 
kr —)• 0 , the functional zs —>• 5rJ7’[Eo](zs') rapidly decays for zs away from zd and has a sharp 
peak when zg —>• with a focal spot size of half a wavelength of the incident wave. Moreover, 
it synthesizes the sensitivity 0 /j 7 ’[Eo](z 5 ) relative to the insertion of an inclusion at zs G ri. 
Heuristically, if the contrast (e 2 r — 1) has the same sign as {eir — 1), then ^[Eo](z 5 ) must 
observe the most pronounced decrease at the potential candidate Z 5 G for the true location zjj. 
In other words, zs ^p^(zs') is expected to attain its most pronounced negative value; refer, 
for instance, to for detailed discussions on sign heuristic. Notice that is positive if 

the contrasts of true and trial inclusions have same signs. Consequently, by virtue of the decay 
property, dTJ\^o]{zs) = — (zg) assumes its maximum positive value when zs —>■ z^. In 
a nutshell, achieves the detection of zd and resolution limit. 


2.1. Imaging with multiple measurement 

Let G be equidistributed directions and 9^’^ {£ = 1,2) be such that ,6^’“^} 

is an orthonormal basis of We define the incident fields and the topological sensitivity 
functional by 

^ 2 n 

E;j’^ := iKOj x "" and dxJizs) ■= - ^ ^ clrJ"[E;^’^](z5). (18) 

^ i=i j=i 


As for n sufficiently large ^ YlJ=i ^ 
therefore 

2 n 


1 JkOJ (x-y) 


: jo(K|x — y|) and ^ basis of 


i ^ X, = i ^ (I _ 

i=\ i=l J=1 


T\ ^ikOJ (x-y) 


i=i 


- ^ (I - ^VxVx ) ~- 47 r(Keo) ^Am{r(x,y)}. 


2 n 


-47r(/^eo) ^9m{r(x,y)}. 


(19) 


Similarly, 

- IL 

£=i i=i 
Following result holds. 

Theorem 2.4. Let zg G n gN be sufficiently large and A : B = Then 

drJizs) « {T{zs, zd)} Mp : {r(z5', z^)} | + 0(/)- (20) 

Proof. Since A9 ■ 9 = A : 99'^, the approximation (|19l) yields 

3 3^ 2 n _ 

^;^^y]y]5Re{9m{r(z5,z^)}MpE^’^(zD) •M5E^’^(z5)| +0(p^), 


drJ (zs) = - 

p^K^C. 


i=i i=i 


dvreon 


EE 5Re{Am{r(zs,z^)}Mp : Ms{9j x 9f’^)i9j x + O(p^), 


£=1 j=l 


p^K'^eQ^CMe{Qm{T{zs,ZD)}Mp : {r(zs,z^)} } +0{p‘^). 


□ 


Theorem 12.41 substantiates that drUizs) oc 9m{r(zs, z/j)} and has a peak (sharper than 
that of c1tx7[Eo]) when zs —)• z^) (see Remark 1 2.3 p . 







3. Stability with respect to measurement noise 

The aim in this section is to establish stability of the imaging functional BtJ- Assume that the 
measurements of the far held arnplitude are corrupted by a mean-zero circular Gaussian noise 
^ —>• C^, that is, E“(x) := E^(x) -|- ^(x) for all x G where E“ is the corrupted value. 

Henceforth, a superposed ~ indicates true value without noise corruption. Further, let ci^ be 
the noise covariance of ^ such that for all y, y' G 


E 


^(y)^(y') = )i 


and 


E 


(y') 




( 21 ) 


where j and j' indicate the j-th and /-th measurements and bjy is Kronecker’s delta. 

The noise affects c1ti7[Eo] and dTJ by means of Herglotz wave, = 'H_e:[E“] +'He[^]- 

The term 'H_e:[E“] is independent of noise and yields true image; discussed earlier. On the other 
hand, T-LeI^] is a circular Gaussian random process with mean-zero by linearity and by dehnition 
of Furthermore, for all z, z' G 


E 




=E 

=E 


§2x§2 


^(x)^(y)'^e 




ds{x)ds{y) 

X X (,t(x) X ds{±)ds{y) 


/§2xS2 

To get the latter identity, the fact that ^ G Ly(S^) is used. By virtue of (fT^ - (fT3]l . one gets 

=E [ 

O JS2x§2 


E 




(I-ii^)^(k)^(y)^e*"" ^ ds(k)ds(y) 


(l- [ e'"*^("“"')ds(k) = -^^9m{r(z,z')}. 

\ ) Js2 


( 22 ) 

Thus, T-LeI^] precipitates a speckle pattern, that is, a random cloud of hot spots having typical 
diameters of the order of wavelength and amplitudes of the order of / ^/k. The covariance of 
corrupted image is furnished by Theorem 12.21 and expression (1221) as 

Cov{dTJ{'Z‘),dTJ{'Z‘')) 

2 2 n _ _ 

=J E E {^E[e%) • M,Ej^(z)} )Re • M5 E^’''(z')} 

2 2 n f _ 

=£2 E E M^Ef (z) • E \HE[e^]{^) [nE[e^]{'^'] 


t=i j=i 






re^KCo 


• 3m{r(z,z')}M5Ef (z')} 


i=i 


where = K^(47r) ^ ^€ 2 ^ — l) and the identities ([2T]) - ([22]) are invoked. 

Now substituting expression for E^’^ and using approximation (I19p . we arrive at 

Cov(9tJ"(z),<9t1^(z')) 

O 2 2 2 n 

= - ^ X ef ’). am{r(z, z')}m,(«, x , 

” i=l j=l 

^K‘^cr^^{e2r - l)^(2n)“^eQ^3fte|M59m{r(z,z')} : 9m{r(z, z')}M5|, 



































In order to fathom the statistics of corrupted image, assume for an instance that is a sphere 
so that M ,5 = 3{2€2r + Then, it can be readily verified that 

Cov{dTJiz),dTJ{z')) « ||9'm{r(z, z')} ||^ , (23) 

where = 3\Os\{t2r ~ + 1) ll-^ll •“ VA : A is the Frobenius norm of A. 

Remark 3.1. The above analysis elucidates that noise perturbation of image is typically of order 
a^lV^ and the hot spots have shape identical with that of actual peak procured by dxJ subject 
to accurate measurements. The main spike of OtJ is not altered. Further, the dependence 
of typical perturbation size on Xj^fTn indicates that the topological sensitivity based imaging 
framework is more stable when multiple measurements are available at hand. 

3.1. Signal-to-noise ratio 

It follows immediately from (j23p that the variance of 8x31 at z G can be approximated by 

Var(9T^(z)) ~ 6^(T^^K"^(2n)“^||ijm{r(z,z)}||^. (24) 

Specifically when Dp is spherical inclusion, the signal-to-noise ratio (SNR) defined by SNR = 
E [dxJi'ZD)] (Var [dTJ{'2‘D)])~^^‘^ , can be approximated by virtue of (1241) and Theorem 12.41 as 

SNR« 3V^a^-^eo^|2er/ + - l\p^\OD\n^pm{T{zD,ZD)}\\. (25) 

The expression (I25p indicates that signal-to-noise ratio depends directly on contrast, volume of 
the inclusion, the operating frequency and inversely on noise covariance. 

4. Conclusion 

The topological sensitivity functions presented herein for detecting a dielectric inclusion, from 
single or multiple measurements of electric far field amplitude, are proved to be very efficient in 
terms of resolution, stability and signal-to-noise ratio. In fact, the Rayleigh’s resolution limit is 
achieved and stability with respect to measurement noise is substantiated. The stability with 
respect to medium noise is not dealt with, however, moderate stability can be guaranteed under 
Born approximation. The reader is referred, for instance, to mm in this regard. 
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